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Abstract 

Mass  transport  phenomena  through  the  membrane  have  the  major  role  in  performance  of  a  proton  exchange  membrane  (PEM)  fuel  cell.  Water  is 
transported  through  the  membrane  due  to  diffusion,  convection  and  electro-osmotic  drag.  Normally  it  is  assumed  that  the  membrane  is  impermeable 
against  gases.  Strictly  speaking  there  are  some  amount  of  solute  gases  in  liquid  water  which  move  within  the  membrane  and  reach  to  the  other 
side.  For  example,  oxygen  dissolves  into  the  water  at  cathode  catalyst.  Most  of  this  oxygen  reacts  with  proton  and  produces  electricity.  But  some 
oxygen  molecules  diffuse  toward  anode  side  and  react  directly  with  hydrogen  at  anode  catalyst.  This  process  leads  to  waste  of  energy  since  direct 
reaction  releases  energy  in  form  of  heat.  In  this  paper  these  processes  have  been  studied.  We  have  developed  a  two-dimensional  numerical  model 
using  full  Navier-Stokes  equations  and  species  transport  equations  of  hydrogen,  oxygen  and  water.  It  was  found  that  cross-over  of  reactant  gases 
has  a  considerable  effect  on  cell  energy  efficiency  in  some  common  cases. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  study  deals  with  many  physical  and  theoretical 
issues.  The  fluid  is  a  mixture,  containing  various  species  with 
completely  different  properties.  The  species  movement  is  gov¬ 
erned  by  several  mechanisms.  Diffusion  and  convection  are  com¬ 
mon  transport  mechanisms  for  all.  For  water,  electro-osmotic 
drag  within  the  membrane  is  very  important.  At  liquid-gas  phase 
interfaces  we  have  dissolution  of  gases  such  as  hydrogen,  oxy¬ 
gen  and  nitrogen. 

Moreover  chemical  reactions,  which  occur  in  a  very  thin  lay¬ 
ers  (catalysts)  together  with  changes  in  porous  media  properties, 
affect  the  flow  pattern,  concentrations  and  hence  fluid  properties 
especially  in  regions  closer  to  catalysts. 

Strong  coupling  between  flow  parameters  is  another  issue. 
For  instance  concentration  filed  of  species  depends  on  fluid 
velocity.  In  the  other  hand,  local  mechanical  properties  of  fluid 
such  as  density  are  a  weighted  function  of  species  properties  and 
hence  depend  on  each  specie  concentration. 
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In  such  complex  problem  analytical  models  can  not  be  so 
popular.  Their  scope  of  work  is  limited  to  give  estimation  in 
very  simplified  cases  neglecting  many  effective  aspects.  It  must 
be  mentioned  that  analytical  methods  have  their  own  rank  since 
they  give  parametric  results. 

In  the  other  hand,  experimental  methods  can  give  reliable 
results.  However,  they  have  drawbacks  such  as  high  costs  and 
measuring  limits.  Many  local  quantities  can  not  be  measured 
inside  the  cell,  because  the  dimensions  are  so  small  and  the  cell 
performance  is  intensely  sensitive  to  sensors. 

Considering  above  mentioned  limitations  of  analytical  and 
experimental  methods,  Numerical  methods  have  received  most 
attention  in  fuel  cell  researches.  These  methods  can  give  various 
information  in  a  wide  variety  of  operating  conditions. 

Many  numerical  works  have  been  accomplished  during  last 
15  years.  One-dimensional  isothermal  model  of  Bernardi  and 
Verbruge  [1]  which  contained  only  cathode  was  a  base-frame 
for  the  next  investigations.  Singh  et  al.  [2]  presented  a  two- 
dimensional  model  considering  cross-sectional  area  normal  to 
membrane  as  the  solution  domain.  They  used  Darcy  model 
instead  of  Navier-  Stocks  and  emphasized  on  better  modelling 
of  mass  transfer.  Kermani  et  al.  [3]  investigated  importance 
of  energy  equation  for  prediction  of  liquid  water  formation. 
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Nomenclature 

Clyj 

water  activity 

c 

concentration  (molem-3) 

Cp 

specific  heat  capacity  (Jkg-1  K-1) 

D 

mass  diffusion  coefficient  (m2  s-1) 

£° 

open  circuit  cell  voltage 

F 

faraday  constant  (C  mol-1) 

GDL 

gas  diffusion  layer 

G 

mass  flux  (kg  m  2  s-1)  or  (mol  m  2  s-1) 

Gcr 

mass  flux  of  gas  crossing  over  through  the  mem¬ 
brane  (mol  m-2  s-1) 

H 

Henry’s  law  constant  (Pam3  mol-1) 

^reac 

liquid  water  enthalpy  of  formation  (J  mol-1) 

i 

current  density  vector  (Am-2) 

icr 

cross-over  equivalent  current  density  (Am-2) 

j 

transfer  current  density  (Am-3) 

k 

thermal  conductivity  (W  m-1  K-1) 

K 

permeability  of  porous  media  (m2) 

Lx ,  Ly 

dimensions  see  Fig.  1  (m) 

M 

molar  weight  (kg  mol-1) 

membrane  molar  weight  (kg  mole-1) 

^eod 

electro-osmotic  drag  coefficient 

P 

pressure  (kPa) 

P  sat 

saturation  pressure  of  water  (kPa) 

R 

gas  constant  (J  mol-1  K-1) 

RH 

relative  humidity 

^O-N 

volumetric  ratio  of  O2  to  N2 

Su 

non-homogeneous  source  term 

Sp 

homogeneous  source  term 

T 

temperature  (K) 

u 

velocity  vector  (ms-1) 

ycell 

actual  cell  voltage  (V) 

yact 

activation  over-potential  (V) 

X 

concentration  (mol  fraction) 

x,y 

coordinates  (m) 

£ 

porosity 

0 

phase  potential  (V) 

A 

water  content  (mol  water/mol  S03_) 

M 

viscosity  (kgm-1  s-1) 

P 

density  (kgm-3) 

Pm 

dry  membrane  density  (kgm-3) 

a 

proton  conductivity  (£2-1  m-1) 

£imP 

impermeability  efficiency 

superscripts 

D 

diffusion 

DR 

direct  reaction 

eod 

electro-osmotic  drag 

g 

gas 

R 

main  reaction 

ref 

reference 

S 

solute,  dissolution 

subscripts 
a  anode 

c  cathode 

eff  effective 

f  fluid 

k  species  (hydrogen,  oxygen,  water  or  nitrogen) 

s  solid  media 


In  recent  years  three-dimensional  models  have  been  developed 
[4-9]  to  investigate  various  aspects  of  processes  occurring  in 
cell. 

Nearly  in  all  works  it  is  assumed  that  membrane  is  an  imper¬ 
meable  wall  against  gas.  It  means  that  only  water  and  proton 
are  present  in  membrane  and  reactant  gases  (hydrogen  and  oxy¬ 
gen)  never  meet  each  other.  So  the  only  occurring  reactions  are 
hydrogen  dissociation  (H2^2H+  +  2e-)  at  anode  and  water 
production  (2H+  +  1/202  +  2e-  — >►  H2O)  at  cathode. 

Impermeability  against  gas  is  an  ideal  property  for  membrane 
as  an  electrolyte,  but  unfortunately  gases  can  dissolve  into  the 
liquid  water  at  GDL-membrane  interface  and  move  to  the  other 
side  of  membrane  [10]  by  diffusion  and  convection.  This  effect 
is  more  important  when  speaking  about  reactant  gases,  because 
they  can  meet  each  other  directly  and  react  at  catalysts.  We  call 
this  phenomenon  as  “secondary  reaction”  or  “direct  reaction”. 
In  fact  direct  reaction  it- self  has  not  a  significant  influence  in 
fuel  cell  performance.  The  notable  issue  could  be  loss  of  fuel 
(and  oxygen  if  pure  oxygen  is  used). 
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Fig.  1.  Solution  domain. 
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Direct  reaction  at  anode  and  cathode  depends  on  amount  of 
oxygen  and  hydrogen  cross-over  respectively.  Reactant  gases 
solubility  in  liquid  water  is  small.  Diffusion  coefficients  of  gases 
in  dissolved  form  are  a  few  orders  of  magnitude  less  than  that 
in  gas  phase.  Also  because  of  greater  value  of  density  and  low 
porosity  of  media  in  membrane  the  velocity  is  negligible.  So  it 
can  be  concluded  that  normally  the  amount  of  gas  cross-over  is 
low  and  major  portion  of  reaction  occurs  in  indirect  form  (main 
reaction)  to  generate  electricity. 

Here  the  importance  of  cross-over  is  investigated  from  point 
of  view  of  reactant  gases  consumption. 

2.  Model  description 

As  Fig.  1  shows,  solution  domain  contains  five  regions: 
gas  diffusion  layer  (GDL)  (anode  and  cathode),  catalyst  layers 
(anode  and  cathode)  and  membrane.  Gas  channels  and  bipolar 
plates  are  considered  as  boundaries. 

The  model  is  developed  under  following  main  assumptions: 

-  The  cell  works  under  steady  state  conditions. 

-  The  fluid  in  GDL  is  ideal  gas.  Volume  of  liquid  water  in  this 
region  is  negligible  and  liquid  phase  has  no  effect  on  gas  flow. 

-  The  flow  all  over  the  cell  is  laminar. 

-  Product  water  is  in  liquid  phase. 

-  All  crossed-over  reactant  gases  are  consumed  in  direct  reac¬ 
tion,  i.e.  neither  hydrogen  exists  in  cathode  GDL  nor  oxygen 
in  anode  GDL. 

-  GDL  has  infinite  resistance  against  proton  flow.  So  all  protons 
generated  at  anode  catalyst  migrate  to  cathode  catalyst  and 
react  with  oxygen. 

-  Porous  media  has  the  same  temperature  as  the  fluid. 


Spd  is  source  term  modelling  friction  between  fluid  and 


porous  media  (Darcy’s  model): 


2 

h  £  U 

Spd  = 


(5) 


y 


It  is  worth  saying  that  the  velocity  used  in  the  equations 
is  only  due  to  convective  mass  transfer  and  diffusive  veloc¬ 
ity  which  originates  from  chaotic  motions  of  molecules  is  not 
included. 


2.1.3.  Energy  conservation 

d<y(>CpT)  +  V.(epCpuT)  =  V.deffVr)  +  Suenergy 
3 1 


(6) 


Energy  equation  is  applied  for  both  porous  media  and  fluid. 
Consequently,  treatment  of  interactions  between  these  two 
phases  is  not  needed,  but  the  properties  must  be  calculated  as 
the  average  value  of  solid  media  and  fluid. 

Source  term  consists  of  ohmic  losses  of  proton  current,  acti¬ 
vation  over-potential  and  direct  reaction: 


/  -2 
r 


- 1-  (LaCtJ  +  G^Y.h react) 

a 

i* 
a 

;2 

ract  •  i  Cr 


r>Cr 


Suenergy  =  J 


a 


anode  —  catalyst 


membrane 


- h  (V  j  +  G  ./z react )c  cathode  —  catalyst 


cr 


0 


otherlocations 


(7) 


i  is  current  density  magnitude. 


2.1.  Governing  equations 

2.1.1.  Continuity 

d2ld  +  v.[s(pu  +  Gd)]  =  Sumass 
3 1 


Gd  is  defined  as  diffusive  mass  flux: 


gd  =  £-  Acvxkpk 


(i) 


(2) 


Su  encompasses  mass  source/sink  terms  due  to  indirect  reac¬ 
tions  and  electro-osmotic  drag  as  is  introduced  for  source  terms 
of  species  transport  equation: 


(Su)mass  =  Y,  f(S<  +  Su|od)pk 


i  concentraion 


(3) 


2.1.2.  Momentum 

3  (spu)  _>  j 

- 1-  V .(spun)  =  —sVP  +  V.(e/xeffVw)  +  Sp  (4) 


2.1.4.  Species  transport 


9(g*k) 
3 1 


+  V.(suXk)  =  V.(DkVXk)  +  Su£ 


concentration 


(8) 


Subscript  k  denotes  one  of  species  (hydrogen,  oxygen,  water 
and  probably  nitrogen).  Generally  the  source  term  includes 
consumption/production  of  species  in  main  reaction  and  direct 
reaction  and  also  specie  transport  under  electro-osmotic  drag 
or  dissolution  effects.  Nature  of  source  term  concerns  to  corre¬ 
sponding  specie  as  follows: 

Main  reaction: 


j 


Sl'H: 


2  FC 


tot 


0 


anode  catalyst 
other  locations 


SA 


—  4F'C 1 1  cathode  catalyst 
0  other  locations 


(9) 


i 


SuR 

^  uWater 


7FC  tot 

0 


cathode  catalyst 


3 1 


other  locations 
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Direct  reaction: 


Sugf 


Su™ 


Cry  DR 

Water 


dG%/dx 


dG^/dx 


c 


tot 


0 


cathode  catalyst 
other  locations 

anode  catalyst 
other  locations 


(  2(dGf)/(dx) 


=  < 


Cto, 

yCr 


— (3Gjjp/3x 


C 


tot 


lo 


anode  catalyst 

cathode  catalyst 
other  locations 


(10) 


Ctot  denotes  total  volumetric  mole  of  existing  fluid.  G^  is 
oxygen  flow  rates  crossing  over  through  membrane  and  con¬ 
tains  both  diffusion  and  convection.  So  the  term  — (3Gq^)/(3v) 

denotes  oxygen  consumption  rate.  Similarly  G^  is  hydrogen 
flow  rate  from  anode  site  to  cathode  catalyst.  Note  that  the  term 
dG^/dx  is  inherently  negative. 

Electro-osmotic  drag: 


Vmeod 

ECtot 


(11) 


i  is  current  density  vector  indicating  both  magnitude  and  direc¬ 
tion  of  protons  migration. 

Dissolving: 

When  a  gas  is  in  contact  with  liquid  some  molecules  penetrate 
in  the  liquid  and  dissolve.  The  rate  of  dissolving  is  proportional 
to  gas  partial  pressure: 


Ax  is  length  of  corresponding  control  volume  in  v  direction. 
This  term  is  present  at  gas-liquid  interface  and  shows  the  rate  of 
matter  leaving  gas  region  and  entering  liquid  so  it  has  a  negative 
value  in  gas  phase  region. 

It  is  evident  that  the  dissolution  source  term  is  applied  only 
for  hydrogen  and  oxygen  (and  nitrogen  if  its  transport  equation 
is  included). 

Species  transport  equation  is  solved  for  three  species  hydro¬ 
gen,  oxygen  and  water.  Nitrogen  concentration  is  then  obtained 
knowing  that  the  summation  of  all  concentrations  is  equal  to 
unity: 


An2  =  1  —  Ar2  —  Xo2  —  A  Water 

(16) 

2.1.5.  Phase  potential  (proton  transport ) 

V(crV0)  =  Su0 

(17) 

where: 

f+j  anode  catalyst 

Su0  =  <  —j  cathode  catalyst 

(18) 

[  0  other  locations 

2.2.  Boundary  conditions 

At  v  =  0  (anode  channel-GDF  interface): 

RH^.P^  | 

XWa,er(0,  y)  =  ‘  ^ 

P( o,  v) 

(19) 

VlVr  V)  =  1  -  Xwater(0,  }’) 

(20) 

Xo2(0,  y)  =  Xn2(0,  y)  =  0 

(21) 

At  x-Lx  (cathode  GDL-channel  interface): 


Gm  =  KxPk  =  KxX\P  (j^) 


(12) 


Pk  denotes  partial  pressure  of  specie  k.  superscript  g  is  for 
emphasis  on  gas  phase. 

Simultaneously  some  dissolved  molecules  pass  through  the 
liquid  surface  and  escape  to  gas  phase.  This  can  be  formulated 
as  follows: 


A\Vater(Tv,  y)  — 


RHr.P< 


sat 


P( o,  y) 


X02  (Tx,  y)  —  (1  Awater(0,  y))  X 


An2(Tv,  y)  =  (1  —  Xwater(0,  j))  X 


^O-N 

1  +  7?o-n 

1 

1  +  Rq-n 


(22) 

(23) 

(24) 


Gout 


=  K2Csk 


(13) 


superscript  s  denotes  the  concentration  of  dissolved  form  and 
CjS  is  the  concentration  of  dissolved  specie  in  liquid  adjacent  to 
gas  interface. 

The  value  of  dissolved  gas  in  liquid  surface  falls  at  equilib¬ 
rium  point  which  the  entering  rate  is  equal  with  exiting  rate: 

KlX{P=K2Cl^Cl=  ^-X{P  (14) 

Where  H-K^IKx  is  Henry’s  law  constant. 

Instead  of  direct  calculation  of  concentration  in  liquid  (above 
equation),  dissolution  can  be  considered  numerically  as  a  source 
term: 


£uS  KlXjP-K2Cl 


(15) 


XU2(Lx,  y)  =  0 


(25) 


At  bipolar  plates  interfaces  we  have  zero  mass  flux. 

For  energy  equation  convective  heat  flux  was  considered  at 
channel  interfaces. 

For  phase  potential  equation,  GDF  was  excluded  from  solu¬ 
tion  domain.  As  the  proton  flux  at  GDF-eatalyst  interface  is  zero 
we  have: 


30 

dx 


GDL — cataly  stint  erfaces 


(26) 


At  top-end  and  bottom-end  of  the  cell,  periodic  boundary 


conditions  were  used  for  all  equations: 

A(x,  o+)  =  A(x,  Ly ) 

(27) 

A(x,  Ly~)  =  A(x,  0) 

(28) 

where  A  is  an  arbitrary  variable. 

GtotAv 
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2.3.  Fluid  properties  and  electrochemical  quantities 


Density: 

P  =  Ctot)  ^xkMk 

k 

Viscosity: 


pkiJk  Y,kx^M^k 
P  ~  Ek^kMk 


effective  viscosity  for  porous  media  is  [9] : 


Meff  — 


1.5(1  -s) 


n  2 


fi 


(29) 


(30) 


(31) 


Thermal  conductivity: 

Effective  thermal  conductivity  is  calculated  as  a  combination 
of  fluid  and  porous  solid  thermal  conductivities: 

1 

Ke  ff  =  -2KS  + -  (32) 

s/(2Ks  +  Kf)  +  (1  -  s)/3Ks 

Kf  =  JfxkKk  (33) 

k 

Heat  capacity: 

Heat  capacity  is  also  considered  as  the  average  value  of  fluid 
and  porous  solid  heat  capacities: 


pCp  =  s(pCp) f  +  (1  -  s)(pCp)s 


(, pCp)f  =  ffpkCpk  =  CxolffxkMkCpk 

k  k 


Mass  diffusion  coefficients: 
For  gas  phase  we  have: 


(34) 

(35) 


(36) 


In  dissolved  form  diffusion  coefficient  obeys  another  relation 

[11]: 


'  f(T10  A  <  2 

lCT10(-3  +  2A)  2  <  A  <  3 

Dt  =  {  ]n 

10-10(6.5  —  1.167A)  3  <  A  <  4.5 

1.25  x  10-10  4.5  <  A 

v  — 

X  is  water  content  in  membrane. 

Proton  conductivity  [4,6]: 

a  =  (0.5139A  -  0.326)c1268«1/303)-(1/r» 


(39) 


(40) 


Water  content  in  membrane  [4,5]: 

{0.043  +  17.18  aw  —  39.85  a ^  +  36.0  a ^  0  <  a  <  1 
14.0  +  1.4(«w  —  1)  1  <  a  <  3 

16.8  3  <  a 

(41) 


Activity  of  water: 


<2W  — 


water 


sat 


Water  concentration  in  membrane: 


C 


water 


=  A 


Pm 


M 


m 


Total  mole  of  matter  per  volume: 
In  gas  phase: 


Ctot  = 


RT 


In  liquid  phase: 


Ctot  = 


E  Water 
A  Water 


Electro-osmotic  drag  coefficient  [5,7]: 
«eod  =  0.0029A2  +  0.05A  -  3.4  x  10~19 


Current  density  in  membrane: 


(42) 


(43) 


(44a) 


(44b) 


(45) 


D*  =  c- 


T 

[i 


(37) 


\i  denotes  solvent  (water)  viscosity,  c  for  each  solute  gas  is  a 
constant. 

For  liquid  water1  [5,7]: 


Dl  =  Bp  2416((l/303)— (1/D) 

^ water  L 


(38) 


1  The  D\  expression  for  3  <  X  <  4.5  has  been  corrected  here.  In  some  references 
(for  example  [5]  and  [7]),  the  values  of  D,  in  this  range  are  in  contradiction  with 
the  original  equation: 


J  3.1  x  10“7A(e°-m  -  \)e~234-6T  0  <  A  <  3 

(  4.17  x  10“8A(1  +  l6ex)e~2346T  3  <  A 


i  =  -dV  0 


Transfer  current  density: 


J  a  — 


•ref 

7a 


Ot2 

y^ref 

CH2 


0.5 


(e 


Kac,af/RT 


a  _ 


T/act^^/RT 

eva  «a  ) 


_  -ref  [  ^Q2 
7C  7  c  1  ^ref 

co2 


(e 


Vjctar/RT 


c  - 


Vactar/RT\ 

e  c  c  ) 


Cell  voltage: 


yCQll  _  g0  _ 


a 


lx 


o 


dx  +  vp  +  vc 


act 


(46) 


(47) 

(48) 


(49) 
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2.4.  Scaling  cross-over 


Table  1 

Properties*1  and  physical  parameters 


The  amount  of  reactant  gases  cross-over  can  be  expressed 
with  an  imaginary  current  density  named  cross-over  equivalent 
current  density,  /cr: 


=  if  +  if 

(50) 

=  4  FG%2 

(51) 

=  2  FGg2 

(52) 

if  shows  rate  of  oxygen,  passing  from  cathode  to  anode  through 
the  membrane.  In  the  other  hand  if  all  crossing  over  oxygen 
reacts  with  proton  at  cathode  catalyst  instead  of  direct  reaction 
at  anode  catalyst,  then  we  will  produce  if  with  this  amount  of 
oxygen.  Similar  statement  can  be  said  for  hydrogen  and  if. 

It  is  useful  to  have  a  quantitative  statement  for  energy  effi¬ 
ciency  of  the  cell  in  terms  of  impermeability  against  reactant 
gases  as  follows: 

For  hydrogen: 


.imp  __ 
ML  - 


cg2 


+  g£2 


For  oxygen: 


+  /ave 


.imp  _ 
*02  “ 


cr 

go2 


+  cg2 


;cr 

a 


if  +  zave 


Overall  impermeability  efficiency: 


(53) 


(54) 


parameter 


Unit 


Value 


s  in  GDL 

s  in  membrane  and  catalyst 

Volume  fraction  of  membrane  in  catalyst 

K  in  GDL  (permeability  to  air) 

K  in  membrane  (hydraulic  permeability) 

Pm 

Mm 

/ef  at  anode 
/ef  at  cathode 


m 


hh2  =  (t)H2at  323  K 

H°2  =  (i)02at  353  K 
H°2  =  (i)02at  3“  K 


2 

m 

kg  m-3 
kg  mol- 
Am-2 
Am-2 


Pa  m3  mol 


0.4 

0.25 

0.5 

1.76  x  10-11 
1.8  x  10-18 
1840 
1.1 

5  x  108 
1  x  102 


jTCf 

K 

353 

pref 

kPa 

101.325 

h  reat 

Jkg-1 

-285830 

R 

Jkg-1  K"1 

8.314 

F 

C  mol- 1 

96487 

Dref 

H? 

2  -1 
m  s 

1.1  x  10-4 

Dref 

2  -1 
m  s 

3.2  x  10-5 

n'U 

^Water 

2  -1 
m  s 

7.35  x  10-5 

jy> 

h2 

at 

353  K 

2  —1 
m  s 

1.06  x  10-8 

Ds 

h2 

at 

323  K 

2  -1 
m  s 

3.40  x  10-9 

A 

at 

353  K 

2  —  1 
m  s 

7.47  x  10-9 

do2 

at 

323  K 

2  -1 
m  s 

2.43  x  10-9 

Hh2 

=  1 

(i)H2at  353  K 

Pa  m3  mol- 1 

3.2886  x  10 

2.7426  x  104 
2.0265  x  104 
1.0387  x  104 


a  Values  for  gases  dissolution  and  diffusivity  in  water  have  been  adopted  or 
calculated  from  [11],  porous  media  properties  are  from  [9]  porosity  in  membrane 
has  been  assumed  according  to  [12]. 


£imP  =  1  - 


GCHr2 


cg2 


+  Gl_ 


/>cr 

02  ^  ^imp^imp 

Gg2  +  Gg2  ~  ^  ^ 


2.5.  Numerical  procedure 


(55) 


A  finite  volume  code  was  developed  to  solve  governing  equa¬ 
tions. 

Power  law  scheme  was  applied  to  approximate  diffu¬ 
sion/convection  terms.  The  SIMPLE  algorithm  with  staggered 
grid  system  was  adopted  to  solve  the  equations. 


Table  2 


Operating  conditions  for  the  base  case 


parameter 

Unit 

Value 

Ro-n  at  cathode  channel  entrance 

mol  mol- 1 

0.266 

Fluid  temperature  at  anode  channel 

K 

353 

Fluid  temperature  at  cathode  channel 

K 

353 

Pressure  at  anode  channel  entrance 

atm 

3 

Pressure  at  cathode  channel  entrance 

atm 

3 

RH  at  anode  channel  entrance 

% 

95 

RH  at  cathode  channel  entrance 

% 

50 

3.  Results  and  discussion 

For  evaluation  of  the  present  model,  several  test  cases  were 
studied.  The  physical  parameters  are  listed  in  Table  1.  A  base 
case  with  operating  conditions  and  dimensions  listed  respec¬ 
tively  in  Tables  2  and  3  was  studied  for  comparison  the  simulated 
results  with  experimental  data  and  also  for  overview  of  species 
transport. 

Fig.  2  compares  calculated  polarization  curve  of  the  base 
case  with  the  experimental  data  reported  by  Wang  et  al.  for  a 
similar  operating  conditions  [12].  A  good  agreement  can  be  seen 
between  the  present  results  and  experimental  data. 

Fig.  3  gives  current  density  contours  of  the  base  case.  This 
figure  illustrates  that  the  current  density  in  x-direction  is  approx¬ 
imately  10  times  greater  than  that  in  y-direction.  It  means  that 
protons  almost  move  straight  toward  the  cathode  catalyst. 


In  Fig.  4  hydrogen  concentration  is  displayed.  The  minimum 
concentration  within  the  anode  GDL  is  located  in  adjacent  with 
catalyst,  closed  to  axis  y  =  Ly  where  we  have  most  distance  from 
channels.  Since  hydrogen  should  be  dissolved  in  water  to  enter 
membrane  region,  a  discontinuity  at  anode  GDL-membrane 
interface  is  occurred.  Fig.  5  shows  oxygen  concentrations.  Simi- 


Table  3 

Dimensions  and  corresponding  grids  for  the  base  case 


Region 

Length  (fxm) 

Number  of  allocated  grids 

Gas  diffusion  layers  thickness 

200 

17 

Catalysts  thickness 

20 

1 

Membrane  thickness 

108 

10 

Cell  width  (Ly) 

2000 

61 

bipolar  plates  width 

800 

— 
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larly  a  discontinuity  is  observed  when  oxygen  goes  to  membrane 
at  cathode  catalyst.  The  minimum  concentration  of  oxygen  in 
cathode  GDL  is  located  at  y  =  Ly ,  next  to  catalyst  boundary. 

With  the  above  background  in  proton  transport  and  reactant 
gases  distribution,  we  can  study  gas  cross-over  (direct  reaction) 
contribution  in  waste  of  energy. 

In  the  first  step  we  consider  the  effect  of  current  density.  For 
this  reason,  several  tests  were  carried  out  with  similar  conditions 
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Fig.  3.  Current  density  within  the  membrane  in  (a)  x-direction  (b)  y-direction. 
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Fig.  4.  Hydrogen  concentration  (mol  fraction). 


of  the  base  case  for  a  wide  range  of  current  densities.  The  results 
are  presented  in  Fig.  6.  According  to  this  figure,  icr  decreases 
as  current  density  increases.  This  is  due  to  more  reactant  gas 
consumptions  at  greater  current  densities.  For  example,  hydro¬ 
gen  concentration  is  decreased  at  anode  side  of  the  membrane 
if  its  consumption  (dissociation  in  anode  catalyst)  is  increased, 
so  fewer  hydrogen  molecules  have  opportunity  to  pass  through 
the  catalyst  and  diffuse  to  cathode. 

In  the  next  step  three  effective  parameters  were  considered 
in  four  cases  as  listed  in  Table  4.  The  average  current  density  in 
all  cases  is  5000  Am-2.  The  calculated  results  are  also  given  in 
this  table. 

Case  I  corresponds  to  the  minimum  gas  cross-over.  Com¬ 
paring  with  case  II  (base  case)  effect  of  membrane  thickness  is 
observed.  Thickness  of  membrane  acts  as  a  resistance  against 
both  gases  and  proton.  Consequently  with  a  thicker  membrane, 
the  rate  of  gas  cross-over  through  membrane  decreases  and 
ohmic  loss  increases.  It  must  be  noted  that  even  though  a  thicker 
membrane  leads  to  more  impermeability  efficiency,  it  has  more 
ohmic  losses  and  this  losses  usually  play  a  greater  role  in  loss 
of  output  power. 

From  case  III  it  is  evident  that  use  of  pure  oxygen  as  oxidant 
instead  of  air  leads  to  a  significant  increase  in  oxygen  cross¬ 
over  and  strengthens  anode  direct  reaction  in  comparison  with 
cathode  side.  It  can  be  seen  that  a  380%  increase  in  oxygen 
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Table  4 


Study  cases  and  corresponding  results 


Case 

Problem  conditions 

Results 

Membrane  thickness  (pm) 

Fluid  in  cathode  channel 

Channels  fluid  temperature  (K) 

/“(Am-2) 

(“(Am-2) 

£imP  (%) 

I 

230 

Air  with  50%  relative  humidity 

353 

5.0 

11.3 

99.7 

II  (base) 

108 

Air  with  50%  relative  humidity 

353 

11.7 

26.2 

99.2 

III 

108 

Pure  oxygen  with  50%  relative  humidity 

353 

64.2 

26.2 

98.2 

IV 

108 

Pure  oxygen  with  50%  relative  humidity 

323 

75.4 

19.7 

98.1 

concentration  from  mole  fraction  of  0.191  (case  II,  standard  air 
with  RH  =  50%)  to  0.915  (case  III,  pure  oxygen  with  RH  =  50%) 
has  led  to  a  450%  increase  in  oxygen  cross-over.  In  this  case  rate 
of  energy  loss  has  almost  amounted  to  the  considerable  value  of 
2%. 

Case  IV  highlights  the  effect  of  temperature.  Temperature 
affects  two  main  parameters  in  gas  transport  through  mem¬ 
brane:  solubility  and  diffusivity.  As  the  solubility  of  gases  in 
liquid  water  is  greater  in  lower  temperatures,  their  diffusiv¬ 
ity  becomes  smaller.  The  resulting  cross-over  can  increase  or 
decrease  depending  on  contribution  of  each  parameter.  In  case 
IV  the  amount  of  hydrogen  cross-over  is  less  than  that  of  case 
III,  but  due  to  greater  oxygen  cross-over,  the  total  losing  current 
density  has  slightly  grown. 

Although  direct  reaction  is  negligible  at  high  current  density 
conditions  for  cells  with  thick  membrane,  we  observed  that  in 
some  cases  it  is  considerable,  especially  at  low  current  density 
conditions  which  are  very  common  in  fuel  cells  operation.  In  fact 
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Fig.  5.  Oxygen  concentration  (mol  fraction). 
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Fig.  6.  Cross-over  equivalent  current  density  and  impermeability  efficiency  ver¬ 
sus  current  density. 

even  if  the  cell  is  at  no-load  conditions,  direct  reaction  won’t  be 
stopped  and  even  is  intensified. 

In  the  other  hand,  when  the  cell  is  working  at  a  low  cur¬ 
rent  density,  it  gives  more  voltage,  so  with  a  distinct  cross-over 
equivalent  current  density,  power  losses  (i.e.  product  of  equiva¬ 
lent  losing  current  and  voltage)  increases. 

4.  Conclusions 

A  comprehensive  study  was  performed  to  investigate  the 
effect  of  gas  cross-over  through  the  membrane  in  waste  of 
energy.  It  was  found  out  that  this  effect  is  not  always  negli¬ 
gible  especially  in  calculation  of  cell  efficiency.  The  following 
items  are  concluded: 

-  As  the  membrane  is  thinner,  direct  reaction  increases. 

-  Direct  reaction  occurs  mainly  at  cathode  catalyst  at  moderate 
operating  temperature  when  air  is  used  as  the  oxidant.  With 
pure  oxygen,  direct  reaction  at  anode  obviously  increases. 

-  Direct  reaction  becomes  more  important  at  low  current  density 
conditions. 
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